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A liquid surface touching a solid usually deforms in a near-wall meniscus region. In this work, we replace part of the free 
surface with a soft polymer and examine the shape of this 'elasto-capillary meniscus', result of the interplay between elasticity, 
capillarity and hydrostatic pressure. We focus particularly on the extraction threshold for the soft object. Indeed, we demonstrate 
both experimentally and theoretically the existence of a limit height of liquid tenable before breakdown of the compound, and 
extraction of the object. Such an extraction force is known since Laplace and Gay-Lussac, but only in the context of rigid 
floating objects. We revisit this classical problem by adding the elastic ingredient and predict the extraction force in terms of 
the strip elastic properties. It is finally shown that the critical force be increased with elasticity, as is commonplace in adhesion 
phenomena. 
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Because interfaces are sticky, it is more difficult to pull out an 
object resting on a liquid surface than to simply lift it. This 
is a common illustration of capillary adhesion. The excess 
force actually corresponds to the weight of the liquid column 
drawn behind the pulled object, but relies intimately on sur- 
face tension as the liquid column is sculpted by capillarity 
into a meniscus. Laplace in 1805 was the first to investigate 
this problem in his founding monograph on capillarity M The 
theory designed by Laplace predicted a maximal height for 
the liquid meniscus above which no more solution could be 
found, hence a maximal force needed to extract an object. On 
Laplace's request Gay-Lussac undertook careful experiments 
on the traction of glass disks floating on water and obtained a 
beautiful agreement with this first ever theory of capillarity. 

Adhesion is not the only consequence of the existence of 
menisci. Collective behaviors such as the clustering of bub- 
bles into rafts are also related to this capillary interface de- 
formation. Such capillary attraction (or repulsion) has been 
known for long, and even used in the early 50 's as a model 
for atomic interactions^. In the case of anisotropic particles, 
clustering becomes directional and allows for self-assembly or 
precise alignment of particles^. 

In all these cases, there is an unavoidable elastic defor- 
mation imparted by the meniscus to the neighboring solid 
substrate, which may become non negligible for soft enough 
materials. Quite recently there has been a keen interest on 
the possibility to strongly deform an object with such elasto- 
capillary interactions. As a matter of fact, the variety of elasto- 
capillary driven shape distortion is surprisingly vast: clump- 
in g0, buckling^, wrinkling^ or selectable packing 8 , as well as 




Fig. 1 Elastocapillary meniscus experiment. A thin polymer 
strip (violet), of length L and width W, is clamped inside a nar- 
row rigid box partially filled with water (blue). Initially, the water 
level matches the embedding height. Upon withdrawal of the liquid, 
the soft strip and the free surface deform altogether into an elasto- 
capillary meniscus of height H. 



changing in wetting properties 9 to quote a few examples, all 
being particularly relevant at small scales^. 

Therefore, in contrast to Laplace and Gay-Lussac experi- 
ment, elastic structures resting on liquid interfaces may ex- 
perience large deformations when dipped in or withdrawn 
from the liquid, because of interactions with hydrostatic pres- 
sure and surface tension. Nature offers an elegant example 
of such deformations with aquatic flowers avoiding flooding 
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Fig. 2 Typical elasto-capillary meniscus shapes, (a) Equilibrium configurations for a Mylar strip floating on water. The length of the strip 
here is L = 21 mm. From left to right, the meniscus height is respectively H = 5.0, 7.2, 8.9 and 12.5 mm. (b) Same snapshot sequence for a 
PVS strip resting on water. The length of the strip is L = 15.8 mm. From left to right, the meniscus height is respectively H = 10.7, 13.6, 16.1 
and 16.6 mm. 



of the pistil by folding their corolla as the level of water in- 
creases 11 . This mechanism recently inspired Reis et al. 12 
a smart technique to grab water from a bath using a passive 
'elasto-pipette'. A different but related problem is the case of 
a thin polymeric film draping over a point of contact on wa- 
ter^, showing the appearance of wrinkles and their transition 
towards localized folds. Elasticity is also the key to understand 
some surprising phenomena involving free surfaces, as for in- 
stance the difficulty to submerge hairs into a liquid by push 
ing thenJSI, the water-repellent ability of striders' legsl 15 l 16 lo: 



or 

even the possibility for some insects and reptiles to walk on 
the water^. On this topic, Burton and Bush^S proposed re- 
cently a simple model to describe how elasticity modifies the 
buoyancy properties of a floating rod. 

In the present study we propose a fully non linear de- 
scription of the shape a soft thin object resting on (and de- 
formed by) a liquid surface. We study particularly the quasi- 
static extraction of such an object from the liquid bath and 
we consider the deformed liquid-solid ensemble as an elasto- 
capillary meniscus - a generalization of the well known cap- 
illary meniscus including elasticity at the interface. In the fol- 
lowing we will detail the different regimes appearing in the 
problem, depending on the relative importance of elasticity, 
gravity and capillarity. Particular attention will be paid to the 
behavior of the contact line between liquid and solid, whose 
slippage can seriously influence the extraction process. We fi- 
nally discuss the extraction process at the light of the particular 
regime considered and show the possibility to increase the ex- 
traction force with a purely elastic mechanism, as is common 
in adhesion phenomena. 



2 An elasto-capillary meniscus 

In our experimental setup, sketched in Fig. [T] a narrow 
box with rigid glass walls is filled with water (density p = 
10 3 kg/m 3 and surface tension 7= 72 mN/m). A thin and 
narrow strip, of length L, width W and thickness e, is clamped 
at one end to the box, and is free at the other end. Strips em- 
ployed in the experiments are made of polyethylene terephtha- 
late (PET, bending stiffness per unit width B = 1.7 • 10 -4 Nm) 
or polyvinyl siloxane (PVS, bending stiffness per unit width 
5=1.8- 1(T 5 Nm). The lateral extent of the box is D > W, 
such that a small gap exits between the strip and the walls in 
order to avoid friction. 

The initial amount of water in the box is set to leave the 
liquid surface at the same level than that of the embedding in 
the wall, 7 = 0. The strip floats on water and is completely flat. 
Then, rather than to extract the strip, water is quasi- statically 
withdrawn from the box, so that the liquid interface lies at 
distance H from the embedding. The liquid-air interface pins 
to the free edge of the strip, causing its deflection. 

Focusing on the (X,Y) plane, we show Fig. [2] several equi- 
librium configurations, obtained for one PET and one PVS 
strip, and for increasing values of the parameter H. The over- 
all deformation of the liquid surface and the elastic strip results 
from a competition between capillarity, elasticity and gravity. 
It shall thus be referred to as an elasto-capillary meniscus in 
the following. 

Looking at the shape of the free liquid surface, it is already 
quite clear that the total height of the compound can be much 
larger than that of a pure liquid meniscus. Note that there is 
no paradox however, as the height of the liquid portion still 



Fig. 3 Elastic strip extraction. Sequence illustrating the failure of the elasto-capillary meniscus, associated with air invasion from the end of 
the strip and violent elastic relaxation. The time intervals between two consecutive images are 168 ms, 78 ms and 27 ms. 



does not exceed 2L gc . The elasticity of the strip is the fac- 
tor allowing to reach such heights, and to lift large amount of 
liquid. 

One can observe that the main difference between the se- 
quence in panel [2ja) and[2jb) is related to the evolution of the 
contact angle between the strip and the liquid. Indeed, in panel 
[2|a) the angle decreases with respect to //, from a value near 
180° in the first image to approximatively 120° in the last im- 
age. Conversely in panel [2jb) the contact angle increases with 
respect to H up to approximatively 230° in the last frame. The 
variety of contact angle values is a consequence of the sharp 
edge at the end of the strip. This contact line pinning avoids 
air or water invasion in the region below or above the strip 
respectively, as already pointed out in Reis et alM^, provided 
that the hysteresis of the contact angle is strong enough as we 
will discuss later. 

The height of the elasto-capillary meniscus can be increased 
up to a point where there is a sudden receding of the contact 
line associated with air invasion and a violent elastic relax- 
ation, as depicted Fig. [3] (see also the corresponding movie in 
Supplementary Material). This is the extraction of the strip 
from the interface. The fact that the extraction is not a con- 
tinuous but a violent process, occurring for a critical height 
depending a priori on all the parameters of the problem, will 
be examined in detail in section [4] 



3 Theoretical description 

In this section we present a 2D theoretical description of the 
shape of the elasto-capillary meniscus. On the one hand, the 
2D approach for the elastic strip is justified by the fact that in 
all experiments L ^> W and L^> e, which allows for Euler- 
Bernoulli description of a beam. On the other hand, one can 
notice that the liquid interface is also deformed in the Z di- 
rection, as suggested by the thick dark shape of the interface 
in Fig. [2] This is the consequence of the contact between the 
water and the glass walls. See Supplementary Information for 
a discussion on the role of this 3D effect. 



3.1 Liquid portion 

Pressure everywhere in the liquid is given by hydrostatic law: 

p(Y)=p a -pg(Y + H) (1) 

where p a is the pressure in the air, p the density of the liquid 
and g gravity acceleration (see Fig. [4] for notations). When- 
ever H ^ 0, the liquid portion such that Y > —H is in de- 
pression with respect to the atmosphere. Pressure p(Y) jumps 
to p a across the curved liquid interface, due to Laplace law 
jk = Ap. We introduce the angle y between the tangent to 
the meniscus and the X-direction (i// is positive if trigonomet- 
ric), which is related to the curvature by K = — ^ = —\j/ f (S), 
S being the arc-length along the meniscus. Laplace law can be 
written: 

YV\S)=pg(Y m (S)+H) (2) 

where Y m is the vertical position of the liquid interface. Using 
the differential relation Y^(S) = sin \ff(S), we can finally write: 

y/'(S) = -^sinyA(S), (3) 

with L gc = (y/pg) 1 / 2 the gravity-capillary length. In the case 
of a meniscus vanishing far from the wall, \f/(S — » °o) = 0, 
equation ^ has an analytical solution 19 , which can be writ- 
ten: 

y/(S) = 4 arctan (tan ^ exp — ] . (4) 

Here Yo is a constant, whose value is in general linked to the 
contact angle at the wall. In our system however, the contact 
angle changes with H, as already displayed in the previous 
section. We will show later the way in which x/Zq can be found. 

3.2 Elastic portion 

The deformation of an inextensible beam can be expressed as 
a function of the angle 6 between the tangent to the beam and 
the X-direction (6 is positive if trigonometric). 

As the beam lies at the interface, it is subjected to air 
pressure on the upper face and water pressure on the lower 
face. Thus, the effective distributed force per unit width is 




Fig. 4 Notations used in the model. This sketch is a 2D view of 
Fig. [T] introducing the notations used in the model. The local an- 
gle between the X-axis and the strip (violet) is denoted 0, the one 
between the X-axis and the free surface (blue) is and the contact 
angle between the liquid and the strip is (p. 



f = — pg (Y + H) n, where S is the arc-length along the beam, 
spanning from to L, and n = — sin 6e x + cos 0e y is the uni- 
tary vector normal to the upper face. Note that we neglect the 
role of the own weight of the polymer strip, as it floats on the 
liquid interface. 

Equilibrium of internal force F(5) = F x (S)e x + F y (S)e y and 
internal moment M(5) = M(S) e z is given by Kirchhoff equa- 
tions, whose projections on the three axes are: 

= -pg(Y(S)+H)smG(S) (5a) 
F;(S) = pg(Y(S)+H)cos6(S) (5b) 
BO" (S) = F x (S) sin (S) - F y (S) cos (S) (5c) 

These equations are written per unit width. We introduced the 
bending constitutive relation B6'(S) = M(S), with B = EI /W 
the bending stiffness per unit width (E is the Young modulus 
and / the second moment of the area of the cross section). The 
position of an infinitesimal element ds is related to the angle 
6 by the relations: 



X'(S) = cos G(S) 

Y'(S) = sin 0(S) . 

Finally, boundary conditions are: 

X(0)=0 ; 7(0) =0 
0(0) =0 ; 0'(L)=O 



(6a) 
(6b) 



(7) 



meaning that the strip is clamped horizontally at the origin and 
curvature vanishes at the end. Two more boundary conditions 
are related to the force applied by the liquid-air interface at 
the tip of the strip: here surface tension pulls the beam with 
a force of magnitude y and an angle y/b with respect to the 
horizontal (see Fig. [4]). Hence: 



It is convenient to write these last boundary conditions as a 
function of the contact angle (p between the solid and the 
liquid rather than y/o. Looking at figure [4] it appears that 
\j/o = 0(L) + (p — n, therefore: 



F X (L) = -y cos (6 (L) + q>) 
F y (L) = -ysin(e(L) + p). 



(9a) 
(9b) 



Kirchhoff equations ([5]) can be written in a more compact form 
(see Supplementary Information for details): 

BB"\S) + 1 0\sf + ye\S) cos (p + pg(Y(S) +H) = (10) 

A characteristic length scale appears naturally in this prob- 
lem when a balance between elastic and hydrostatic forces 
applies, the elasto-hydrostatic length L e h = (B/pg) 1 / 4 . This 
length scale was already introduced by Hertz^ and Foppl^H 
in their founding works on elastic plates. The elasto- 
hydrostatic length is also the relevant length scale when con- 
sidering the buckling of an elastic beam resting on a compliant 
substrate^, in which the rigidity of the foundation is replaced 
by pg- We use L e h to make lengths dimensionless: 

{ W ,M,4c} = 7^ {S,XJ,H,L,L gc } (11) 

^eh 

Dividing equation ( [TO] ) by pgL e k leads to the system of 
equilibrium equations (dependence on s in now suppressed for 
brevity): 

Q m + 1 0'3 +£ 2 c CQS 9 0' + (y + A ) = o (12a) 



(12b) 
(12c) 

(13) 



x = cos 

y' = sin 

with boundary conditions: 

jc(0)=0 ; y(0)=0 ; 0(0) =0 
e 7 (£)=0 ; G"(£)=f gc sm(p 

3.3 Matching condition 

Because of contact line pinning at the edge, the contact an- 
gle (p between the elastic strip and the liquid surface does not 
result from classical Young's construction. Rather, (p is set 
by the geometrical constraint expressing the anchoring of the 
meniscus at the strip edge: 



Jstrip — ^meniscus 

(£). 



(14) 



F X (L) = y cos wo ; F y (L) = y sin y/Q . 



(8) 



This sole condition does not guarantee a constant value for <p 
and indeed, the contact angle exhibits experimentally a strong 
dependence with the non-dimensional depth h. 
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Fig. 5 Comparison between experimental and theoretical pro- 
files, (a) The theoretical equilibrium shape (obtained with £ = 1.83, 
£ gc = 0.24 and h = 1.09) is superimposed to its experimental coun- 
terpart (corresponding to the last frame of Fig. [2]panel a), (b) Same 
comparison with the second frame of Fig. [2] panel b (here £ = 2.44, 
£ gc = 0.48 and h = 2.06). 



After integrating equation ^ once, the right hand side of this 
condition can be rewritten as 



^meniscus M = ~2£ gc sin — - h , (15) 

leading to the following convenient end point condition for <p: 

(16) 



); strip (£) = 2£ gc cos( ) _ h 



4 Elasto-capillary meniscus equilibrium 
shapes and their limit of existence 

4.1 Predicted and observed equilibrium shapes 

Typical equilibrium configurations observed in experiments 
for PET and PVS strips are reported Fig. [5] These two exam- 
ples correspond to different values of the three dimensionless 
parameters (h, £ and £ gc ) governing the overall shape of the 
structure. In order to test the assumptions made in the model, 
we confront these observations with theoretical predictions. 
Specifically, we integrate numerically the equations ( [T2|[T3] ) 
associated with the matching condition ( [T6| ). The shape pre- 
dicted with the model for the experimental governing parame- 
ters is then superimposed to the observed profiles, and a good 
agreement can be noted both for PVS and PET strips. 

Is the observed deflection of the strip gravity-driven or 
capillarity-driven? To answer this question we introduce the 
analogue of the elasto-hydrostatic length for capillarity, i.e. 
the elasto-capillary length 5 L ec = (B/y) 1 ^ 2 which is the typ- 
ical length needed for capillary bending effects to overcome 
the flexural rigidity of the strip. As L QC L gc = L^ h we see that 
the nondimensional parameter £ gc can be written as the ratio 
Lqw/Lqc- Hence £ gc is a measure of the relative importance of 
hydrostatic forces over capillarity in the bending process. In 
our experiments, £ gc is systematically lower than 0.5 but still 
of order unity, meaning that hydrostatic pressure is hardly the 
dominant effect in the deformation of the strip. Therefore we 
chose not to neglect £ gc in the following discussion. 
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Fig. 6 Limit of existence for the elasto-capillary meniscus. The 

paths followed by equilibrium states are shown in the (h, <p) plane, 
in which (p is restrained to physical values between and 2n. The 
curves correspond to increasing values of the nondimensional strip 
length £, starting with a near zero value for the black curve (from 
bottom to top £ = 0.01, 1, 1.5, 1.8, 2, 2.2, 2.4, 2.6, 2.8, 2.9, 3, 3.2 
respectively). It appears that for each curve a maximal value for h 
exists, and beyond this limit point no more solution can be found. 
The dashed blue line reports the path of these limit points as £ is 
varied. Three equilibrium shapes, corresponding to the limit points 
Pi , P3 and P5 , are reported on the right. Note that throughout this 
diagram, pinning of the contact line on the strip tip is assumed. 



4.2 The path towards extraction 

We now investigate the nature of the extraction process per 
se, with a particular focus on the dependence of the critical 
height with the governing parameters. To this end, we make 
use of a continuation algorithm to follow the equilibrium con- 
figuration as the strip is withdrawn from the bath. These equi- 
libria branches are reported Fig. [6] Note that for comparison 
purposes the value of the nondimensional capillary length has 
been set throughout this study to £ gc = 0.24, corresponding to 
the experiment involving a strip made of PET. 
We start by considering the rigid limit case £ = 0.01 where 
the elasto-hydrostatic length L e h is 100 times larger than the 
length of the strip L (downmost black curve in Fig[6j. When 
h = 0, the liquid surface is perfectly horizontal = 0) and the 
apparent contact angle is (p = 71. As h is increased, the con- 
tact angle drops down to zero for h = 2£ gc = 0.48. In this rigid 
case, h simply measures the true height of the liquid meniscus. 
This height cannot exceed the critical height 2£ gc correspond- 
ing to a purely wetting state (<p = 0, see point P\) - otherwise 



the matching condition (14) would not be fulfilled. This is 
precisely the critical height obtained by Laplace 1 that allowed 
him to predict the critical extraction force, measured experi- 



mentally by Gay-Lussac. 

Now adding some elasticity in the problem, we consider 
analogous equilibrium branches for flexible strips ranging 
from £ = 0.01 up to £ = 3.2. At first glance, it can be no- 
ticed than whenever h = the deflection of the strip is nil and 
so the contact angle (p is always equal to 71. Starting from this 
flat configuration, we follow each branch of equilibrium as h 
varies. It is noteworthy that the addition of elasticity allows 
to reach higher maximal heights in the system. This is clearly 
illustrated by the elasto-capillary menisci shapes at points P3 
and P$ in Fig. [6] 

Most of the equilibrium branches in Fig. [6] exhibit a limit 
(fold) point. In other words, two equilibria coexist for a given 
value of h\ presumably one stable and the other unstable. At 
the critical value, these two equilibria coalesce and then dis- 
appear: there is no more equilibrium solution past the limit 
point. This behaviour is commonplace in dynamical systems 
(saddle-node or fold bifurcation), and was first noticed in the 
context of interfaces by Plateau 23 . More precisely, when a 
catenoid made of soap film is formed between two rigid rings, 
there is also a limit point where two catenoids (one unsta- 
ble and the other stable) merge and then cease to exist. Ex- 
perimentally, beyond the limit point, the catenoid bursts into 
droplets, and the elasto-capillary meniscus breaks down {i.e. 
the lamella is extracted) as illustrated Fig. [3] Note that for 
the two uppermost curves (£ = 3 and £ = 3.2) the maximum 
value for h no more corresponds to a limit point, but rather to 
a constrained maximum occurring when (p = 271, the largest 
possible value for the contact angle. 

The locus of these maxima, be them limit fold points or 
constrained maxima, is revealed Fig. [6] with a dashed blue 
line. The path followed by these maxima appears non-trivial, 
suggesting optimum values for the elasticity to achieve largest 
heights and exhibiting a tumble as (p approches 271. As these 
points mark the frontier where equilibria disappear, it is there- 
fore quite natural to expect from this line to represent the ex- 
traction limit for the strip. 

4.3 Pinning and slipping of the contact line 

A simple way to verify the exactitude of the extraction predic- 
tion is to compare the theoretical shape of the elasto-capillary 
meniscus to its corresponding experimental counterpart. 

Fig. [7] (left) we show two experimental pictures taken at 
the moment of the extraction. In the first case (£ = 2.2- top 
left image), the experimental contact angle at extraction accu- 
rately matches the angle 94 expected from the theory (94 cor- 
responds to the angle at point P4 in Fig.[6j. This means that the 
folding point P4 accurately captures the threshold for extrac- 
tion, comforting our previous analysis. However, the second 
case shown Fig.[7](£ =1.0- bottom left) exhibits a poor agree- 
ment with the theoretical prediction. There, the experimental 




Fig. 7 Slipping of the contact line. Top left: experimental shape 
observed just before air invasion sets in. The strip has a nondimen- 
sional length £ = 2.2, whose limit point in Fig.[6]is the point P4. The 
superimposed angle is the contact angle reached at the limit point P4, 
and it is in good agreement with the one observed in the experiment. 
Bottom left: snapshot corresponding to a strip with £ = 1. Here, air 
invasion is triggered for a contact angle (p clearly above $2, the con- 
tact angle associated with point P2. The instability involves a slipping 
of the contact line as soon as the angle reaches some critical value (p c . 
Right: Equilibrium state for £ = 3.68. The strip literally plunges in 
the bath, with an effective 'interfacial length' £* < £. 



critical contact angle is much larger than the value (p2 expected 
from the location of folding point P2 in Fig. [6] 

The reason for this discrepancy lies in the contact line pin- 
ning hypothesis, that proves to be faulty for some values of 
the contact angle. Indeed, experimentally we clearly observe 
an advance of the contact line as soon as (p approaches a criti- 
cal value (p c ~ 80°. The slippage of the contact line, associated 
with air invasion underneath the elastic strip, is an unstoppable 
process; it proceeds until the lamella is fully extracted. 

Analogously, when the free surface overturns the strip tip a 
similar advance of the contact line can be observed. Actually 
it is triggered as soon as the air wedge angle is less than <p c , 
that is, when (p > 2n — (p c . But this time, the slippage is not 
necessarily a runaway process leading to extraction, as exem- 
plified Fig. [7] (right) with a stable configuration. There, water 
invasion sets in as soon as the critical angle 2% — (p c is reached. 
But the system now settles to a new 'plunged state' that re- 
duces the effective 'interfacial length' of the strip £*(<£). In 
this type of configuration, part of the strip remains stably un- 
derwater. 

These experimental observations tend to depict a somewhat 
more complex picture of the extraction process than previ- 
ously expected. Specifically, the locus of the folding points 
in Fig. [6] does not necessarily describe the extraction thresh- 
old. Extraction cannot be grasped without paying attention to 
the slipping motions of the contact line on the upper or lower 
side of the strip. 

Contact line slippage can however be included in the model. 
Indeed, assuming that the contact angle is constant during slip- 




Fig. 8 Equilibrium branches including pin-slip effects. Loci of 
equilibria in the (h,(p,(£ — £*)/£) space are shown for £=1.5 (blue 
line) and £ = 3.4 (orange line). As long as the contact angle (p lies 
in the range ((p c ^7t — q> c ) (i.e. between the blue stripes), the con- 
tact line is pinned to the strip tip and equilibrium branches lie in the 
(h, <p) plane. When the contact angle reaches a critical value a bifur- 
cation occurs. There, another branch corresponding to slipping con- 
tact line takes off out-of-plane, as the interfacial length £* no more 
corresponds to the actual length of the strip £. Reconnection to pinned 
states appears to be possible. 



ping^ (and equal to (p c or 2n — <p c ), we can still search for 
equilibrium solutions provided that the interfacial length £* is 
turned into a free parameter. Note that the total number of un- 
knowns in this setting is unchanged as (p is known. The end 



condition (16) can now be used as a condition for £*, and the 



system of leading equations ( [T2|[T3] ) be integrated from to 



Figure [8] displays the equilibrium branches for two values 
of the nondimensional length, i = 1.5 and £ — 3.4, when 
both pinning and slipping of the contact line are retained. 
These branches are plotted in the three-dimensional space 
(A, <p, (£ - £*)/£). Here (£ - £*)/£ represents the relative im- 
mersed/emerged portion of the strip (this portion can either be 
in air or water). Note that this quantity cancels out when the 
contact line is pinned to the strip edge. 

For £ = 1.5 (blue line), the path followed by equilibrium 
solutions is identical to the one already shown Fig. [6] for pure 
pinning. However, as the contact angle (p reaches the critical 
value (p c an out-of-plane bifurcation occurs. The new branch 
arises from contact line slipping along the underside of the 
strip, as shown in the inset. This bifurcation is subcritical and 
the equilibrium branch monotonically tends to (£ — £*)/£ = 1, 
meaning full extraction of the strip (for clarity reasons, the 
branch is only displayed for values such that (£ — £*)/£< 0.3). 
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Fig. 9 Critical meniscus height versus strip length. Comparison 
between the critical height h max reached in experiments as a function 
of the nondimensional strip length £ (orange dots) and the theoretical 
prediction based on the 'pin- slip' end condition (magenta curve). For 
reference, the results using a pure pinning end condition are shown 
with the dashed blue curve. The different regimes in terms of slipping 
of the contact line are also displayed. 



A careful investigation of the nature of the bifurcations occur- 
ring at (p = (p c suggests that those are all subcritical, i.e. lead 
to a sudden extraction of the strip. 

For £ = 3.4 (orange line) the contact angle now increases 
with h. The free surface overturns the strip edge, and a slip- 
ping bifurcation occurs when (p = 271 — (p c . The nature of the 
bifurcation is now different, as equilibrium 'plunged states' 
exist for values of h larger than that at the bifurcation. Typical 
shapes of plunged states are shown in the insets of Fig. [8] and 
illustrate the slipping of the contact line above the strip. The 
loop shape of the out-of-plane branch implies that the inter- 
facial length t first decreases with h until it reaches a mini- 
mum value, and increases up to £ again as shown in the insets. 
In contrast to air invasion the immersion of the strip is thus 
a reversible phenomenon. Surprisingly enough, the slipped 
states branch then exhibits a reconnection with a new in-plane 
(pinned states) branch. It is remarkable to point out that the 
two in-plane orange branches are by no way connected, even 
when (p is allowed to adopt values far beyond its physical lim- 
its (data not shown). This is the ambivalent 'pin-slip' nature 
of the end condition that allows for such a plentiful bifurcation 
diagram allowing to track equilibria along multiple branches, 
at variance with the catenoid case where the end condition is 
unique. 

We now confront the experimental critical height h mSLX 
reached by the elasto-capillary meniscus as a function of £ 
with the model, considering both pure pinning and 'pin-slip' 
conditions (Fig. [9]) (in the pure pinning case the dashed curve 



reporting the critical height is basically a translation of the 
limit points locus shown Fig. [6] in a function of £). Clearly 
the agreement between the experimental observations and the 
'pin- slip' prediction is excellent. The differences between 
pure pinning and pin- slip are even more pronounced for larger 
values of £. Indeed, for pure pinning a strange tumbling be- 
haviour is particularly noticeable in the region 3 < £ < 5, sug- 
gesting the odd idea that lengthening the strip shortens the 
critical height. This wrong prediction is actually linked to 
the non-monotonic shape of the limit points locus. The full 
model incorporating mixed 'pin-slip' end conditions allows 
for a correction of this prediction and smoothes out the evo- 
lution of h max with £. This highlights the importance of con- 
sidering out-of-plane bifurcations and reconnections inherent 
to contact line slipping along the upper side of the strip in the 
extraction process. Note also that Fig. [9] exhibits the whole 
regime phenomenology appearing in this problem, along with 
their non-trivial apparition sequence. Finally, looking at Fig. [9] 
suggests an asymptotic behavior for £ ^> 1 , a possibility we in- 
vestigate next. 

4.4 Universal behavior of long strips 

We now present a theoretical focus on the case £ ^> 1 where the 
strip is very long compared to the elasto-hydrostatic length. 



Rescaling ( 1 2c ) and ( 1 2a ) with s = s/£ and y = y/h leads to: 



£ 



y (s) = sin 0(5) 



(17a) 



^(e-(^) + ie / (,~) 3 ) + |^e / (,~)cos(p = o. (i?b) 



Taking l/£ as a small parameter, we immediately obtain at 
leading order the external solution 6=0 and y = — 1 . This 
external solution represents a strip lying flat over the liquid 
surface, with a contact angle (p = n. A boundary layer solution 
exists and reconciles the external flat solution with the bound- 
ary condition y(0) = 0. We introduce classically the boundary 
layer variable § as follows: s = £§, with § = 0(1) within the 
boundary layer. In terms of this variable, equation ( |17a| ) reads: 



h 1 

£e 



y'(S) = smO(S), 



(18) 



and dominant balance requirements immediately give the 
thickness of the boundary layer as £ ~ h/L 
Injecting this results into equation yields: 



1 



e'"(S)- 



t 1 



(19) 



In this equation the hydrostatic pressure term is of order 
one, and has to be balanced by - at least - one other term. 



Let's first consider the case £ gc <C h, where hydrostatic pres- 
sure is solely balanced by elasticity. The balance can only hap- 
pen if h ~ 1, consistently with the significance of the elasto- 
hydrostatic length L e h. This constraint certainly gives clues 
on the asymptotic bounded behaviour of /z max when £ achieves 
large values, as shown Fig. [9] In this limit (amounting to ne- 
glect capillary forces) the shape of the strip is governed by the 
following equation - expressed with the original variables: 



e /// + ^e /3 +v+/z = o. 



(20) 



A second option in the dominant-balance analysis is to re- 
tain only hydrostatic pressure effects and surface tension ones. 
This corresponds typically to h ~ £ gc and £ gc ^> 1 . The shape 
adopted by the strip in this regime is described by solutions 
of: 

y + h-P gc 0' = 0. (21) 

When derived this equation can be recast into the well 
known meniscus equation 0" = l/i^ c sin0. The shape of the 
elastic strip is thus that of a pure capillary meniscus. Actually 
this is nearly the case, as solutions of this equation cannot meet 



the clamped boundary conditions at the wall, as equation (21 



has lost two derivation orders with respect to equation ( [19] ). A 
third layer, in which elasticity is restored, can however be in- 
troduced to fulfill the correct boundary conditions of the prob- 
lem. 

Eventually, an appealing feature of the full equation ( [T9] ) 
should be pointed out: its structure bears a striking similarity 
with the equilibrium equation of a buckled elastic beam sit- 
ting on a liquid foundation, for which an analytical solution 
exists 25 . We have recently demonstrated that the exact solu- 
tion of equation ( [T9] ) can be obtained as well and is generally 
not symmetrical 



5 Extraction force 

We are now in a position to address the main question raised 
in this study: How does elasticity affect the extraction pro- 
cess? We start by introducing the (dimensionless) pulling 
force / p , that is the vertical component of the force exerted 
by the clamp: 

pi 

f P = -f y (0) = J (y(s)+h) cos d(s)ds -£ 2 gc smxi/ (22) 

Although the pulling force mainly balances the weight of the 
volume of fluid displaced by the strip, it does rely intimately 
on elasticity and capillarity as these effects shape the liquid 
volume. As the imposed height h is increased, the lifted liquid 
volume varies, and so does the pulling force / p . The extrac- 
tion force / ext is therefore defined as the largest value of / p 
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Fig. 10 Extraction force versus elasto-hydrostatic length. Rela- 
tion between the dimensional extraction force per unit width F ext (in 
N/m) and the rigidity of a fixed-length strip, here measured in terms 
of the elasto-hydrostatic length L e h (mm). The length of the strip 
taken here is 27 mm and the gravito-capillary length L gc has been 
set to 2.7 mm. The figure clearly displays an adhesion peak for a 
specific value of L e h around 0.5 L. Representative snapshots of the 
elasto-capillary meniscus at the moment of extraction are disposed 
along the curve. These snapshots correspond to different regimes; 
From left to right: capillary limit, elasto-hydrostatic regime, optimal 
adhesion, rigid limit. In the inset, a rudimentary model where the 
strip is approximated with a linearly varying curvature rod and the 
meniscus with a vertical liquid wall proves to capture the main fea- 
tures of the main graph, both qualitatively and quantitatively. 



with respect to h\ / ext = max^ / p (/i) . Note that in general this 
maximal value is not achieved when h = h max . 

In order to investigate the exact role of elasticity in the ex- 



traction process, we represent Fig.[10]the variation of the force 
needed to pull out a fixed-length strip with variable rigidity. 
More precisely the dimensional extraction force per unit width 
^ext = Pg^eh/ext is represented as a function of L e h for a given 
strip length L. For plotting purposes, we chose a strip length of 
L = 27 mm and a fixed gravito-capillary length L gc = 2.7 mm 
corresponding to our experiments, but the features displayed 



Fig. [TO] appear to be generic. 

Different regimes in the extraction process can readily be 



noticed on Fig. 10 inspection. First, for very small values of 
L e h (corresponding to L e h < L gc , that is £ gc > 1), the strip is 
so soft that it deforms almost with no resistance under the 
combined action of capillarity and gravity, as would have a 
pure liquid interface. This capillary limit is emphasized in the 
leftmost thumbnail of Fig. [H)| with a comparison between the 
elasto-capillary meniscus shape and the looping shape of a real 
liquid meniscus. This excellent agreement is consistent with 
the fact that in the limit £ gc > 1, the elasto-capillary meniscus 



shape is captured by equation (21 



When L e h reaches slightly larger values (but still such that 
L e h <C L, meaning £ ^> 1), elasticity comes into play with the 
elasto-hydrostatic regime. There, the strip is flat and aligned 
with the liquid surface outside a region of size G(HjL^) 
where the deformation is concentrated, as already discussed. 
The shape of the elasto-capillary meniscus in the inner region 
is given by solving equation ( [20} . Note that as long as £ is 
still larger than unity, the shape of the strip exhibits a self- 
similar behavior as it only depends on the elasto-hydrostatic 
length L e h, and so does the lifted liquid volume: F QXi ~ L 2 h 
/ext constant). Stiffening the elastic strip therefore results 
in a quadratic increase of the extraction force. 

For higher values of L e h an adhesion peak clearly enters 
into the picture. This peak, centered at L e h ~ 0.5L (i.e. £ ~ 2), 
corresponds to a cut-off of the self-similar regime imposed by 
end-effects, as L e h becomes of the same order of magnitude as 
L. Here, the large deflection of the strip cause large amount of 
liquid to be drawn behind, hence a maximal extraction force. 
For even larger values of L e h, the elastic strip progressively 
rigidities as the extraction force falls down to its asymptotic 
Laplace/Gay-Lussac limit. It should be pointed out that the 
extraction work, relevant in the context of imposed displace- 
ment, displays a similar adhesion peak (data not shown). 

This sequence of different regimes can be captured with 
a rough model containing as only ingredients elasticity and 
gravity (inset of Fig. [10]). In this model, the curvature of 
the strip is approximated with a linear function ramping from 
H/Ll h at the wall down to at the tip. The liquid meniscus 
touching the strip tip is in turn crudely approximated with a 
vertical liquid wall (see snapshot in the inset). Finding the 
extraction force amounts to a maximization problem: For a 
given L e h, we look for the value of H that maximizes the liq- 
uid volume risen by the strip, taking care that the height of 
the liquid wall cannot exceed 2L gc . This simple model, dis- 
regarding both surface tension forces and nonlinear geometric 
displacements of the strip, recovers nicely the main features of 
the relation F QXt vs the adhesion peak, the different limits, 
the optimal shape. The quantitative agreement of this elemen- 
tary model, free of adjustable parameters, certainly reveals the 
key physical ingredients necessary to understand the role of 
elasticity in capillary adhesion. 

6 Conclusion 

Elasticity is known to increase energetic costs in adhesion phe- 
nomena, and to induce hysteric cycles with different deforma- 
tions back and forth. For example, the long polymer chains in 
elastomers can store all the more elastic energy that the chain 
is long. This energy is lost when the chains rupture and that 
cracks initiate 27 . Analogously, soap bubbles gently put into 
contact can be deformed and remain connected even when the 
imposed separation exceeds the initial critical adhesion dis- 



tanceGHl, thus implying bubbles conformations unexplored pre- 
viously. 

In this study, we have focused on a paradigm setup coupling 
elasticity and capillary adhesion. The elastocapillary menis- 
cus discussed here exhibits all the typical features of elastic 
adhesives: extraction force/work possibly increased by elas- 
ticity, hysteric cycles with deformations of the strip depend- 
ing on whether it touches the interface or not, elastic energy 
stored released violently at extraction. This last observation, 
documented in Fig. [3] and in supplementary material, is also 
reminiscent of cavitation pockets appearing in the traction of 
adhesive films^H 

A key point of our analysis is to highlight the existence 
of optimal strip stiffness corresponding to a maximal extrac- 
tion force (conversely, this force can be as low as the ma- 
terial is soft). This could be exploited for example to im- 
prove the design of soft micro-manipulators or MEMS such 
as elastopipettes 12 . Though the maximal volume displaced by 
the strip does not correspond to the actual amount of liquid 
grabbed from an interface, the extraction force and work do 
have an important role in this problem as they set the limit 
points of extraction. 

Finally, the agreement between theoretical and experimen- 
tal shapes of an elasto-capillary meniscus, along with their de- 
pendence with the elasto-hydrostatic length L e h, suggests new 
metrology tools for a precise measure of the elastic properties 
of soft films. 
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